The relationship between perturbation theory and 
direct calculations of rare-earth transition intensities. 



Michael F. Reid 
Department of Physics and Astronomy 
University of Canterbury, Christchurch, New Zealand 
Gary W. Burdick* 
Department of Physics, Southern College 
Collegedale, TN 37315, U.S.A. 
H. J. Kooy 

Department of Physics, The University of Hong Kong 
Pokfulam Road, Hong Kong. 

September 21 1993 

Abstract 

We use a simplified calculation to demonstrate the equivalence between three 
different methods for calculating transition line strengths. These calculations 
demonstrate the complex interplay between spin-orbit and correlation contribu- 
tions to two-photon transitions in rare-earth ions. 
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1 Introduction 

In a previous paper we investigated many-body perturbation theory calculations of 
two-photon transition intensities in Gd^^ We found incompatibilities between the 
many-body formalism and the earlier calculations of Judd and Pooler and Downer 
and co-workers because the latter contain "unlinked" diagrams. Removal of these 
diagrams seemed to destroy the previous good agreement between theory and exper- 
iment. Agreement was restored, however, when "folded" third-order correlation dia- 
grams involving the Coulomb interaction within the 4/'' configuration were included. 
The question that remained unanswered was how these two apparently incompatible 
calculations could yield similar answers. 

More recently we performed "direct" calculations for Eu^+, in which the actual 
eigenstates of the 4/'' and 4f^5d configurations were used [Q. That work clarified 
the interplay between Coulomb and spin-orbit contributions to two-photon transition 
intensities and provided a way to reconcile the different methods of calculation. 

It is the purpose of this paper to explicitly demonstrate this reconciliation. A key 
point is that the "folded" diagrams represent changes in energy from the configuration 
average. Thus, if we are only interested in a restricted range of energy levels, we can 
use the "true" energies instead, in which case the Judd-Pooler calculations give correct 
results. To show this, we will briefly review the Rayleigh-Schrodinger perturbation the- 
ory. We then consider a simplified four-state system, showing the equivalences between 
direct, many-body perturbation, and Judd-Pooler type perturbation calculations. 



2 Perturbation theory 

The many-body perturbation-theory approach used in Ref. []I| is developed in detail in 
Chapters 9 and 13 of the book by Lindgren and Morrison 0. When the formalism is 
extended to include atom-photon interactions, the energy denominators are modified 
by the inclusion of photon energies 0, but the same diagrams occur. The calculations 
employ the concept of a "model space" . Rather than solving for the eigenvalues and 
eigenstates of the full Hamiltonian, an "effective Hamiltonian" is constructed and di- 
agonalized within the model space, and the expectation values of "effective transition 
operators" are evaluated between the model-space eigenvectors. 

Initially, we use Rayleigh-Schrodinger perturbation theory on the level of Chapter 
9 of Lindgren and Morrison. If we only consider excited configurations of opposite 
parity to the ground configuration, the relevant second- and third-order equations for 
the two-photon absorption transition moments are 

^ {K\r\S){S\r\I) 

^ {K\r\T){T\V\S){S\r\I) 

j^{Ej-ET + oo){Ej-Es + ioy ^ ' 

{K\v\S){S\r\J){J\V\I) 
j^{Ej-Es + u){Ei-Es + uj)- ^ ' 

In these expressions, /, J, and K are states within the model space (4/^) and S 
and T are states in an excited configuration {e.g. 4/^~^5(i) outside the model space. 
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The perturbation, V, includes the spin-orbit and Coulomb interactions, r is the dipole 
moment operator, and u is the photon energy {h = 1). If desired, Eq. (|), may be 
Hermitized by averaging with its conjugate. However, it should be noted that model- 
space calculations are inherently non-Hermitian. 

When excited configurations having the same parity as the ground configuration are 
included, we have the two additional equations, 

{K\r\S){S\r\M){M\V\I) 
fj^,{Ej-Es + u){Ej-EM)' ^ ' 

{K\V\M){M\AS){SW) 
^ {Ei - Em + 2u){Ei -Es + uj)' ^ ' 

where the M are states in an excited configuration of the same parity as the ground 
configuration {e.g. 4f^~^5f). For the simplified calculations that follow, however, only 
a single excited configuration will be considered, and thus contributions from Eqs. 
will not arise. 

The above equations will be used for each of the three calculations which follow. 
The fundamental differences between these calculations are not in the equations used, 
but rather in their differing definitions of the states \i) and energies Ei. For the direct 
calculation, eigenstates and eigenvalues of the complete Hamiltonian will be used, al- 
lowing only contributions from Eq. (|I|). In contrast, the many-body calculation uses 
zero-order eigenstates, necessitating the use of all relevant equations. The Judd-Pooler 
type perturbation calculation uses the zero-order eigenstates within the excited config- 
uration (like the many-body calculation) and uses exact eigenstates within the ground 
configuration (like the direct calculation). Thus, Judd-Pooler type calculations omit 
Eq. (^, which has the perturbation V acting within the ground configuration. In the 
next section, we will show that when careful attention is paid to these differences, all 
three calculations result in equivalent answers. 



3 Simplified calculation 

We consider a simplified system with four states, |a), |c), and \d). The model space 
consists of I a) and while |c) and \d) are excited states. The Hamiltonian, H, is 
divided into a zero-order part, iJo? and a perturbation, V , 



H = Ho + V 






V 







A 

A y 



/ Ca 



V 











Z2 






Z2 



(6) 



where A is the zero-order energy difference between the excited and model spaces 
(the Hartree-Fock energy). The Cj represent the non-spherical part of the Coulomb 
interaction, and the Zj represent the spin-orbit interaction. The effective Hamiltonian 
for the model space is 



H, 



eff 



Zi 



(7) 
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In the general case there could be matrix elements of V {e.g. the odd-parity crystal 
field) connecting \a) and \b) to \c) and \d). In that case, Eq. (|^ would be a first-order 
approximation to H^jf. 

If the Zj are small compared to the Cj, we can take the mixing to first order in 
the Zj and ignore the energy shift. This yields the following expressions for the 
eigenstates of H: 



W) 


= !«) + 1^)77^, 

— 


Ea' 


\b') 


= \b) + \a)-^ATr, 

— 


Ey 


W) 


= \c) + \d)-^ATr, 


Ec' 


\d') 


= M) + |c)t7^, 
— 


Ed' 



Ca , 
Cb , 



The states \a') and \b') are also the eigenstates of H^ff. 
We take the one-photon interaction matrix to be 



(9) 
(10) 

(11) 
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(12) 



where the Tk represent dipole moments. Two-photon transitions between \a') and \h') 
are "spin- forbidden" , as they are allowed only due to the presence of the "spin-orbit" 
matrix elements, Zj, in the Hamiltonian. 



3.1 Direct calculation 

In order to perform an "exact" or "direct" calculation of the two-photon transition 
moment for the "spin-forbidden" \a') to \b') transition, we use exact eigenstates of if, 
rather than those of Hq. This makes the third-order terms superfiuous, and only the 
second-order matrix elements of Eq. (|I|) must be evaluated. 

Rather than write out all possible terms, we present only those that are quadratic 
in T2 and linear in Zi. From Eq. (|l|), we obtain 



T2T2 \ f Z\ 



c„-(A + Crf) + cu; \Ca-Cb 



(13) 



where the first fraction represents the second-order time dependent perturbation, and 
the second fraction gives the proportion of |6) in |a'). Expanding Eq. (|T^) in powers of 
l/(— A + C(j), we obtain 

T2T2Z1 ( 1 , Cd-Ca ^ ^^^^ 



Ca-Cb \-A+LU {-A + iu] 
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3.2 Many-body perturbation calculation 

The states in Eqs. d^-j^) represent configurational states, whereas many-body pertur- 
bation theory is expressed in terms of orbitals. However, the correspondence between 
Eqs. and orbital expressions is exact (see Chapter 13 of 0), and for the pur- 

poses of this discussion, we use Eqs. (|l]"0) to demonstrate the results of a many-body 
calculation. Our calculation is based on Rayleigh-Schrodinger perturbation theory, so 
eigenvalues of Hq are used in the energy denominators. Furthermore, the states in Eqs. 
(|l|-^ are the zero-order states (|i), not In order to obtain the terms containing 

T2T2Z1, we place these "effective transition operators" between the eigenstates of H^ff. 
This leads to the factors Z^jiCa — Cft) in some of the following expressions. From Eq. 
(|I|) we obtain 

T2T2 / Zi \ 

^at^ • ^^^^ 

Restricting V to the spin-orbit interaction in Eq. (0), we only get terms involving Z^- 
However, letting V contain the Coulomb interaction yields 



T2CdT2 ( Zx 



(_A + ^)(-A + ^) \Ca-C, 

Eq. (I) yields 



(16) 



T2T2Z1 T2T2Z1 f Ca — Cb ^ 1^ _^ 

X 7^ Tt). (17) 



(-A + cu)(-A + a;) (-A + c^)(-A + c^) VC„ - C, 

TaTaCb f Z^ \ 

" (-A + cu)(-A + cu) \Ca-Cj ' ^^^^ 

from the spin-orbit and Coulomb contributions, respectively. Adding these terms to- 
gether, we get 

a-al-A-Fa;^ (-A + ^)2 )■ ^ ^ 

which is the same as the "direct" calculation. Note the cancellation between Eqs. (p!7|) 
and ([TsD , which arise from "spin-orbit" and "Coulomb" interactions, respectively. 



3.3 Judd-Pooler perturbation calculation 

The original Judd-Pooler calculation does not include Eq. (H) and thus Eqs. (|1^|18D 
do not arise. However, the denominators now contain eigenvalues of -ffe//, and Eq. (^) 
yields 

\ ( ] (20) 

Coulomb contributions within the excited configuration were considered to be small in 
the calculations of Judd and Pooler, and were therefore neglected. For completeness, 
however, we include them here. Thus, from Eq. (|]) we obtain 



T2CdT2 \ f Z\ 



Xa-A + o^)(a-A + u;); \Ca-C, 



(21) 
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Adding these terms together and expanding in powers of 1/ (—A + u) we obtain 

T2T2Z1 /I Cd-Ca \ 

a^l^AT^ + (-A + c.)2+---j' ^^^^ 

which is the same as the results of the direct calculation (Eq. (|T4|)) and the many-body 
perturbation calculation (Eq. (0)). We have emphasized terms containing T2T2Z1, in 
order to display expressions arising from Eq. (^. In general, when terms involving any 
combination of matrix elements are evaluated, with perturbation expansions carried 
out to n-th order, agreement is exact through the — 1 order of l/(— A + u). Thus, to 
achieve agreement between these three calculations in the third order of l/(— A + u), 
it would be necessary to consider fourth-order perturbation terms. 

In essence, as explained in Ref. the Judd- Pooler calculation has moved the part 
of V that acts within the model space into Hq. This means that the matrix element 
{J\V\I) of Eq. (0) is equal to zero, since there is no longer any V separate from Hq 
acting within the model space. The Judd-Pooler calculation is therefore correct, but 
it must be realized that different energy denominators should be used for transitions 
between different energy levels. 

In Judd-Pooler-type calculations (Eq. ^), terms containing the Coulomb matrix 
elements, Ci, arise from using eigenvalues of Heff, which includes effects of the Coulomb 
interaction. In many-body perturbation calculations, by contrast, they arise only from 
perturbation expressions containing the Coulomb interaction explicitly. Therefore, the 
classification of contributions as "spin-orbit" or "correlation" depends, in part, upon 
the method of calculation. 



4 Many-body perturbation theory for Gd^+ 

For the Gd'^+ calculation discussed in Ref. [|l| the states in Eqs. (|1]-§I) represent states 
in the 4/^ and 4f^5d configurations. Many-body perturbation theory, however, is ex- 
pressed in terms of orbitals. The transformation from Eqs. (|l]-§]) to expressions corre- 
sponding to the diagrams evaluated in Ref. is explained by Lindgren and Morrison 
(Chapter 13, see also 

Ignoring diagrams involving core excitations, we have only one second-order diagram 
and four third-order diagrams involving the spin-orbit interaction, Vso, and the dipole 
moment, r. These are shown in Fig. |l](a-c, e-f), and have corresponding algebraic 
expressions 

4^.E '*'^''"t''' . (23) 

s — fcg -f- 0; 

alail^T- — — T7- — — (24) 



St 



{Ei - et + uj){ei - Es + uj) 
{k\r\s){s\r\j){j\Vso\t) 



E j::r:r_r'.v (25) 



{Ei - es + uj){ej - Es + uj 



ai.ai2_^- -, (26) 

rns {^i ~ + ~ ^m) 

t >r- {k\Vso\m) (m|r|g)(g|r|i) 

O'hO'i > 7 TT T- (27) 

(ei-e^ + 2cj)(ei-e, + cu) 
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In these expressions, i, j, and k label valence (4/) orbitals in our 4/'' model space, s 
and t label virtual orbitals outside the model space {5d) , and the e are orbital energies. 
Each term contains a product of a creation operator and an annihilation operator for 
4/ orbitals. The diagram corresponding to Eq. ( P^D (Fig. |T|(c)) is commonly referred 
to as "folded" §. 

Eq. (|23|) arises directly from Eq. (|l]). Eq. (l2^) arises from Eq. @, where the spin- 
orbit interaction and the dipole moments act on the same electron. "Unlinked" ex- 
pressions (Fig. 11(d)), where the dipole moments and the spin-orbit operator act on 
different electrons, arise from both Eq. (0) and Eq. (^ and cancel. This leaves us with 
the "folded" diagram (Fig. |l](c), Eq. (|25|) ) which arises from Eq. (|^) when the dipole 
moments and spin-orbit operator act on the same electron. 

It is important not to confuse "folded" diagrams with "core excitations" . The folding 
is merely a device to preserve the rules for calculating the denominators [^] (especially 
pp. 304-305). In Fig. |l](c), for example, an electron is excited from state i to j to s 
to k, and the diagram does not imply that there is initially an electron in state j, as is 
clear from the operator ordering of Eq. (pSf ) . 

Eqs. (p6|) and ( ^7|) arise from Eqs. (§) and (^). These have been overlooked in 
earlier work [0, |, 0]. Their possible importance has been recently demonstrated by 
Smentek-Mielczarek ^j. 

Some third-order "correlation" diagrams involving the Coulomb interaction are 
shown in Fig. ^. The folded diagram. Fig. H(e), arises from Eq. (|]), when V is the 
Coulomb interaction. It gives a large contribution in the calculations of Ref . , since 
it represents perturbations to the energy-levels of the configuration caused by the 
Coulomb interaction. An analogous term in our simplified calculation is Eq. (|TB]). As 
in our simplified calculation, it is possible to take these effects into account implicitly 
by using the exact energy levels. Thus, the agreement between the calculation of Ref. 

and the calculation of Judd and Pooler (in which the correlation effects are ab- 
sorbed into the denominators) is not merely fortuitous, but is a result of doing the same 
calculation in different ways. 



5 Conclusion 

Careful application of three different calculation methods (direct [^, many-body per- 
turbation and Judd- Pooler || type perturbation calculations) yields identical 
results. Thus, each of these three methods may legitimately be used in the calculation 
of transition line strengths. 

The Judd-Pooler calculation, while valid, is not compatible with Lindgren and 
Morrison's formulation of many-body perturbation theory. This is because the Judd- 
Pooler zero-order Hamiltonian, Hq, is different when acting upon ground configuration 
states than when acting upon states of the excited configuration, thus destroying the 
cancellation of unlinked diagrams necessary for the many-body perturbation theory. 

In contrast, the direct calculations used in Ref. |Q and illustrated here are compat- 
ible with MBPT, as long as we define the model space to include both the 4/^ and 
the 4/^~^5(i configurations. In that case it is possible to add contributions from other 
excited configurations or from other potential terms, using the techniques of MBPT. 
In fact, we feel that this will be a most profitable direction for future calculations, par- 
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ticularly for examination of transition intensities for divalent lanthanides, where both 
the Judd-Pooler formahsm and the more straight-forward MBPT calculations, using 
4/^ as the model space, have been shown to be insufficient IQ. 

Acknowledgements 

We thank B.R. Judd and L. Smentek-Mielczarek for helpful discussions. 

References 

G. W. Burdick and M. F. Reid, Phys. Rev. Lett. 70 (1993) 2491. 
B. R. Judd and D. R. Pooler, J. Phys. C 15 (1982) 591. 



[2] 
[3] 

[4] 

[5] 

[6] 
[7] 



M. C. Downer, in Laser Spectroscopy of Solids 11, edited by W. M. Yen, page 29, 
Springer- Verlag, Berlin, 1989. 

G. W. Burdick, H. J. Kooy, and M. F. Reid, J. Phys.: Condens. Matter 5 (1993) 
L323. 

I. Lindgren and J. Morrison, Atomic Many-Body Theory, Springer- Verlag, Berlin, 
1982. 

M. S. Pindzola and H. P. Kelly, Phys. Rev. A 11 (1975) 1543. 
B. Ng and D. J. Newman, J. Chem. Phys. 83 (1985) 1758. 
L. Smentek-Mielczarek, Phys. Rev. Lett. (1993), submitted. 



February 1, 2008 



9 



Figure 1: Second and third-order one-electron two-photon absorption diagrams. For 
details of notation see Ref. [^]. A single up-going arrow represents a virtual orbital, a 
double arrow a valence orbital. Folded valence lines have a circle around their arrows. 
Photons are represented by wavy lines and the spin-orbit interaction by a triangle. 
Diagrams (a), (b), (c), (e), and (f) correspond to Eqs. (H), (|2|), (H), (||), and ^ 
of the text, respectively. Diagram (c) is folded, diagram (d) is unlinked. 
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Figure 2: Some third-order two-electron two-photon absorption diagrams. Notation is 
the same as Fig. |1|. The Coulomb interaction is represented by a dotted line. Folded 
diagram (e) dominates the correlation calculation of Ref . [|l| . 
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